Divergence of the Classical Trajectories and Weak Localization 
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We study the weak localization correction (WLC) to transport coefficients of a system of electrons 
in a static long-range potential (e.g. an antidot array or ballistic cavity). We found that the weak 
localization correction to the current response is delayed by the large time Ie = lnh\, where A 
is the Lyapunov exponent. In the semiclassical regime tE is much larger than the transport lifetime. 
Thus, the fundamental characteristic of the classical chaotic motion, Lyapunov exponent, may be 
found by measuring the frequency or temperature dependence of WLC. 

PACS numbers; 73.20.Fz,03.65.Sq, 05.45.+b 



I. INTRODUCTION 



An electron system in a static potential is character- 
ized by the following linear scales: the geometrical size of 
the system, L; the transport mean free path Itr = vpTtr 
being the characteristic distance at which a particle can 
travel before the direction of its momentum is random- 
ized; the characteristic scale the potential energy changes 
over, a; and de Broglie wavelength \p, (for the Fermi sys- 
tem Xp = fi/pp, with pp = mvp being the Fermi momen- 
tum). In the most important metallic regime Xp <C L,ltr- 
The scale of the potential a may be arbitrary and depend- 
ing upon this scale two regimes can be distinguished: 

i) Quantum chaos (QC), > Xpltr- 

ii) Quantum disorder (QD), a? < Xpltr- 

The physics behind this distinction is quite transparent: 
after an electron interacts with the scatterer of the size 
a, the quantum uncertainty in the direction of its mo- 
mentum 60 is of the order of 66 ~ Xp /a. Therefore, the 
uncertainty in the position of the particle 6x on the next. 
scatterer can be estimated as 6x ~ ltr69 ~ ItrXp/aSi 
If 6x <C a, the quantum uncertainty in the position of 
the particle is not important and its motion can be de- 
scribed by the classical Hamilton (or Liouville) equations. 
Except some special cases, these equations are not inte- 
grable, the electron trajectory is extremely sensitive to 
the initial conditions and the classical motion is chaotic. 
The quantum phenomena in such regime still bear es- 
sential features of the classical motion; it is accepted in 
the literature to call such regime "quantum chaos" . In 
the opposite limit, 6x ^ a and the electron looses any 
memory about its classical trajectory already after the 
first scattering. Any disordered system where the Born 
approximation is applicable may serve as an example of 
QD regime. 

Under assumption of the ergodicity of the system, the 
classical correlator is usually found from the Boltzmann 
or diffusion equations. The form of these equations is 
identical for both regimes. The only difference appears in 
the expression for the cross-section entering into the colli- 
sion integral. For the QC, this cross-section can be found 



by solving the classical equations of motion, whereas in 
the QD it is determined by solution of the corresponding 
quantum mechanical scattering problem. 

Subject of weak localization (WL) theory is the study 
of the first order in Xp/ltr corrections to the transport 
coefficients of the system. The WL in the quantum 
disordec, kave been studied for more than fifteen years 
alreadyo The regimg of. the quantum chaos attracted 
attention only recentlydllj. This interest was motivated 
mostly by technological advances which allowed the fab- 
rication of the structures where a 3> A^. Two examples 
of these structures are: (1) the antidot arraysO where role 
of a is pJgyed by the diameter of an antidot; (2) ballistic 
cavitieaj'Q where a kr — L coincides with the size of 
the cavity. 

Weak localization corrections are known to have 
anomalous dependence upon the frequency u, temper- 
ature or the applied magnetic field and that is why they 
can be experimentally observed. For the two-dimensional 
system case L — > oo the WL correction to the conductiv- 
ity Act can be conveniently written as 



r(tj)ln 



1 



UJTtr 



ijJTtr ^ 1, 



(1.1) 



where s = 2 is the spin degeneracy, and r(w) is a renor- 
malization function. It is this function where the differ- 
ence between the quantum disorder and quantum chaos is 
drastic. Gorkov, Larkin and KhmelnitskiicI showed that, 
for the whole frequency domain, F — 1 for the quan- 
tum disorder and does not depend upon the details of 
the scattering. The question is: Does such a universality 
persist for the quantum chaos too? 

In this paper, we will show that, in the limit cj — > 0, 
the renormalization function F — > 1 which proves the 
universalitji-|Of weak localization correction for the quan- 
tum chaostll. However, unlike for the quantum disorder, 
F acquires the frequency dependence at oj much smaller 
than l/Ttr- This frequency dependence can be related 
to the Lyapunov exponent A characterizing the classical 
motion of the particle. It gives an opportunity to extract 
the value of the Lyapunov exponent from the measure- 
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ments of the frequency dependence of the conductivity. 
We found 



r(w) = exp ( 2iujtE 



2uj^X2tE 

A2 



(1.2) 



where Ehrenfest time Ie is the time it takes for the mini- 
mal wave packet to spread over the distance of the order 
of a and it is given byOJ 



1, fa 



(1.3) 



Quantity A2 — A in Eq. (1.2) characterizes the deviation 



of the Lyapunov exponents, and it will be explained in 
Sec. m in more details. In the time representation, result 
( |1.2D corresponds to the delay of the weak localization 
correction to the current response by large time 2tE, see 
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FIG. 1. The weak localization correction to the con- 
ductivity in the time domain, Aa{t) = J ^Aa{uj)e~^"* 
for the quantum chaos (solid line) and quantum disor- 
der (dashed line) regimes. The developed theory is valid 
for t > tE- Solid curve is calculated for parameters 
A = 4A2 = l/rt„ ln(a/AF) = 7. 

The paper is arranged as follows. In Sec . O, we present 
the phenomenological derivation of Eq. ( |L2D . The ex- 
plicit expression relating the weak localization correction 
to the solution of the Liouville equation will be derived 
in Sec. In Sec. IV, we will find the quantum correc- 
tions to the conductivity in the infinite chaotic system. 
Sec. describes the effects of the magnetic field and fi- 
nite phase relaxation time on the renormalization func- 
tion. The conductance of the ballistic cavities is studied 



in Sec. VI. Our findings are summarized in Conclusion. 



II. QUALITATIVE DISCUSSION 

Classical diffusion equation is based on the assumption 
that at long time scales an electron looses any memory 
about its previous experience. However, during its travel, 
the electron may traverse the same spatial region and be 



affected by the same scatterer more than once. These two 
scattering events are usually considered independently, 
because with the dominant probability the electron en- 
ters this region having completely different momentum. 

However, if we wish to find the probability Wq{T,po) 
for a particle to have the momentum opposite to the ini- 
tial one, p(T) — — p(0), (time T is much larger than 
Ttr) and to approach its starting point at small distance 
|r(T) — r(0)| — po a, we should take into account the 
fact that the motion of the particle at the initial and 
final stages are correlated. This is because the trajec- 
tory along which the particle moves on the final stage, 
[r(T — t), p(T — t)] almost coincides with the trajectory 
particle moved along at the initial stage, [r(t),p(t)], see 
Fig. 1^. These correlations break down the description of 
this problem by the diffusion equation. The behavior of 
the distribution function for this case can be related to 
the Lyapunov exponent, and we now turn to the discus- 
sion of such a relation. (Relevance of Wo(T,po) to the 
weak localization correction will become clear shortly.) 




FIG. 2. The classical trajectory corresponding to the prob- 
ability of return at the initial point with the momentum op- 
posite to the initial one. In the "Lyapunov region" the initial 
"i-1" and final "2-f" fragments of the trajectory are governed 
by the same potential. 

The correlation of the motion at the final and initial 
stages can be conveniently characterized by two functions 

pit)=vit)-riT-t), kit)=p{T-t)+p{t). (2.1) 

The classical equations of motion for these functions are 



dp _ k{t) 

dt m ' 
dk _ dU[r{T - t)] _ dU[r{t)] 
dt dv dv 



(2.2a) 
(2.2b) 
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where U is the potential energy. If the distance p is much 
larger than t he c haracteristic spatial scale of the poten- 
tial a, Eqs. (2.2) lead to the usual result {p{t)) cx t^/^ 
at times t much larger than Ttr- Situation is different, 
however, for p <ti a, where the diffusion equation is not 
applicable (we will call this region of the phase space the 
"Lyapunov region"). Thus, the calculation of function 
Wo(T,po) should be performed in two steps. First, we 
have to calculate the conditional probability W{a, po;t), 
which is defined so that the probability for the distance 
p{t) to become larger than a during the time interval 
[t, t + At] is equal to W{a, po;t)At under the condition 
p(0) = po- Second, we have to obtain the probability 
Woia^t) for the diffusively moving particle to approach 
its starting point to the distance of the order of a (it cor- 
responds to the fragment "1-3-2" in Fig. Then, the 
function Wo{T,pq) is given by 



Wo(r,po) 



dtWoia^T -2t)W{a,po;t). (2.3) 



Now, we perform the first step: finding of the prob- 
ability W{a, po',t). We consider more general quantity 
W(p , po; t ) for p < a. We expand the right hand side of 
Eq. (2.2b) up to the first order in p which yields 



dkj_ 
dt 



-M.,j{t)p„ Mi.j{t) 



d''U[v{t)] 
dridri 



(2.4) 



It is easily seen from Eq. (|]^) and Fig. that the change 
in the momentum k during the scattering event is propor- 
tional to the distance p. On the other hand, it follows 



from Eq. (2.2a) that the change in the value of p be- 
tween scattering events is proportional to k. Therefore, 
one can expect that the distance p grows exponentially 
with time. In Appendix ^ we explicitly solve the model 
of weak dilute scatterers Itr S> a and find the expression 
for the distribution function W{p) = {S{t — t{p))), where 
(...) means average over directions of p. Here we present 
qualitative arguments which enable us to establish the 
form of the function W for the general case. 

We notice that, if matrix Ai{t) does not depend on 
time, the solution of Eqs. (2.2a) and (2.4) is readily avail- 
able: 



Pit) 0^ p(0)e 



At 



(2.5) 



where the quantity A is related to the maximal negative 
eigenvalue of A^. We will loosely call A the Lyapun ov ex - 
ponent. If Ai varies with time, the solution of Eqs. (2.2a) 
and (2.4) is not possible. We argue, however, that for the 
large time t Ttr, this variation may be described by a 
random correction to the Lyapunov exponent: 



dlnp 
dt 



= X + SX{t). 



(2.6) 



At time scale larger than Ttr the correlation between the 
values of SX{t) at different moments of time can be ne- 
glected, {SX{ti)6X{t2)) = X2S{ti — t2), that immediately 
gives the log-normal form for the function W: 



W{p,po;t) = 



A3 



27rA2£(p) 



exp 



X{Cip)~Xtf 
2X2C{p) 

In p/pq. 



(2.7) 



Formula ( p.Tl) is valid in general case even though ana- 
lytic calculation of the values of A and A2 (as well as of 
the diffusion constant) can be performed only for some 
special cases, e.g. for Itr ^ a. For the antidot arrays, 
A is given by the inverse scattering time up to the fac- 
tor of the order of ln{ltr/a)£3 The model of the dilute 
weak scatterers is considered in Appendix The result 
is A, A2 — Ttr^ [Itr / aY^^ . In the ballistic billiards, coef- 
ficients A, A2 are of the order of the inverse flying time 
across the system. 



Equation (2.7) describes the distribution function only 
in the vicinity of its maximum, \\n [p / po) — Xt\ 5, Xt. 



However, this result will be sufficient if time T in Eq. (2.3) 
is large enough T > C{a)/2X. At smaller times the prob- 
ability of return is determined by the tail of the distribu- 
tion function W{p) which is by no means log-normal. 

It is worth mentioning, that there is some arbitrari- 
ness in our choice of the initial conditions p(T) = — p(0), 
|r(T) - r(0)| = Po. The other possible choice is |p(r) 



p(0)| = fco, |r(T) - r(0)| = 0. In this case, formula ( p/7| ) 
remains valid upon the substitution po — > ako/p{0). 

Now we can find Wo{T,po) from Eq. dJ). Substitut- 
ing Eq. (2.7) into Eq. ( |2.3D , we arrive to the result for 
the probability Wo (T, po): 



VFo(T,po) 



duj 
2^ 



Wo{uj,po)e' 



'iujT 



(2.8) 



\ ( ^ /2iW'C(a) 2c^2^2/:(a) 
Wo(w,po) = M^D(w,a)exp ( ^3 



where Wu{^tCl) is the Fourier transform of the func- 
tion WD{t,a). Function M^£)(a,a;) = VF£)(ti;;a 
hr)WD{ijJ,ltr) is determined by two consecutive pro- 
cesses. First process, with the probability H^£)(w;a 
Itr), is the separation of the trajectories from distance a, 
at which they become independent to the distance larger 
than Itr, where the diffusion equation is applicable. The 
characteristic time for such process is of the order of Ttr, 
and thus W£)(w;a Itr) = 1 + 0{uJTtr). The probabil- 
ity W£)(w, Itr) is found by solving the standard diffusion 
equation. For the two-dimensional case, which will be 
most interesting for us, function WD{uj,a) has the form 



WDiuj,a) 



1 



■In 



/ 1 



AttD \ LUTtr 



(2.9) 



where D = VpTtr/2 is the diffusion constant. Notice that 
this function does not depend on a. Expressions (2.9) 
and (2.7) are written with the logarithmic accuracy. 

So far, we considered a purely classical problem. We 
found the probability for a particle, propagating in a clas- 
sical disordered potential, to approach its starting point 
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with the momentum opposite to its initial one. In the 
calculation of the classical kinetic coefficients (e.g. con- 
ductivity), the integration over all the direction of the 
momentum is performed. As the result, the peculiarities 
in the probability discussed above are washed out and 
do not appear in the classical kinetic coefficients. How- 
ever, the function Wo{p,t) plays very important role in 
the semiclassical approach to some quantum mechanical 
problems._pne of such problems arose long time ago in 
the studyllJ of break down of the method of the quasi,- 
classical trajectories in the superconductivity theoryEj. 
Another problem is the weak localization in the quan- 
tum chaos and we turn to the study of this phenomenon 
now. 



a) 



It is well knownHOEIl that the probability w for the 
particle to get from, say, point i to point /, see Fig. ^a, 
can be obtained by, first, finding the quasiclassical ampli- 
tudes Aa for different paths connecting the points, and, 
then, by squaring the modulus of their sum: 



= J2\A^f + J2A^A; 



(2.10) 



The first term in Eq. ( 2.10 ) is nothing else but the sum of 
the classical probabilities of the different paths, and the 
second term is due to the quantum mechanical interfer- 
ence of the different amplitudes. For generic pairs a, 
the product AaA*^ oscillates strongly on the scale of the 
order of Xp as the function of the position of point /. 
This is because the lengths of the paths a and /3 are sub- 
stantially different. Because all the measurable quantities 
are averaged on the scale much larger than Xp, such os- 
cillating contributions can be neglected. There are pairs 
of paths, however, which are coherent. The example of 
such paths is shown in Fig. ||b, (paths 1 and 2). These 
paths almost always coincide. The only difference is that 
fragment BEB is traversed in the opposite directions by 
trajectories 1 and 2. In the absence of the magnetic field 
and the spin-orbit interaction, the phases of the ampli- 
tudes Ai and A2 are equal because the lengths of the 
trajectories are close. The region, where the distance 
between trajectories 1,2 is largest, see inset in Fig. 
deserves some discussion. At this point the directions of 
the paths at points Bi , B2 are almost opposite to those 
at points B[,B'2. Furthermore, the differences between 
lengths of paths 1 and 2 should not be larger than Xp- It 
imposes certain restriction on angle Scj) at which trajec- 
tory 2 can intersect itself and on distance Sp to which the 
trajectory 1 can approach itself. Simple geometric con- 
sideration, self-evident from inset i n Fig. p| b, gives the 
estimate Sp ~ V^rhr and Scj) ~ y/Xp/ltr, so that the 
uncertainty relation S(f>dp ~ Af holds. In other words, 
one of the trajectories should almost "graze itself" at the 
point B. 





FIG. 3. Examples of the classical (a) non-coherent and (b) 
coherent paths between points i and /. The scatterers are 
not shown and the paths are straighten for clarity. Encircled 
is the Lyapunov region. The region of the quantum switch 
between trajectories (marked by the rectangular) is blown up 
on the inset. 

The interference part of the contributio n of the coher- 
ent pairs to the probability w, see Eq. ( 2.10 ) is of the 
same order as the classical probability for these trajecto- 
ries. Therefore, the contribution of the interference effect 
to the conductivity a is proportional to the probability to 
find the trajectories similar to those from Fig. ^d. In or- 
der to calculate this probability, we use function Wq {p, t) 
defined in the beginning of this section: the probability 
dP for a trajectory to graze itself during the time interval 
[ti,ti -f dti] is 



dPi = 5p5(t>VpdtiWQ{yJXpltr,ti) = 

XpvpdtiWo{\/ Xpltr, h). 



(2.11) 



in two dimensions. We are, however, interested in the 
correction to the transport coefficients (such as the dif- 
fusion constant or the conductivity). These quantities 
are contributed mostly by the points i, / located at the 
distance ~ Itr from each other. Thus, in order to con- 
tribute to the diffusion constant or the conductivity, ends 
of the trajectories should separate from each other to the 
distance of the order of a, i.e. the trajectories should 
overcome the Lyapunov region one more time. The con- 
ditional probability dP2 that the trajectories diverge at 
the distance ~ a during the time interval [i2j*2 + '^^2] 
under the condition that the self grazing occurred at mo- 
ment ii is given by 
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dP2 = dt2W{a, ^Xphr, h - h 



(2.12) ir,(a;;ri,r2,r3,r4) = Gf+.(ri,r2)Gf_.(r3,r4). (3.1) 



where W is given by Eq. (2.7). 

Summing over all the time intervals, we obtain for the 
quantum correction to the conductivity Act: 



a 



(2.13) 



VfXfJ dt2 J dtiWo(^\/XFkr,tl)^ W(a,^/Xpltr,t2^ . 

If the correction at fi nite frequency uj is needed, the time 
integration in Eq. (2.13) should be replaced with the 
Fourier transform over the total time of travel between 
points initial and final points t = 2t2 + ti in Eq. (2.13). 
This yields 



(2.14) 



where v is t he de nsity of states per one s pin. Coeffi- 
cient in Eq. ( |2.14| ) and sign in Eqs. (|2.13|) and ( ^.14[ ), 
known for the quantum disorder, will be reproduce d fo r 
the quantum chaos in Sec. [II. Substituting Eqs. (2.7) 
and (2.9) into Eq. (2.14) and using the Einst eiH relation 
a = se^D, we arrive to the final result ( |l.2| ).E^ 



III. WEAK LOCALIZATION IN THE QUANTUM 
CHAOS 



It follows from the previous discussion that the calcu- 
lation of the quantum correction is related to the proba- 
bility to find a classical trajectory with largftjisoM'elated 
segments. Standard diagrammatic techniqueuQila is not 
convenient for this case because the averaging over the 
disorder potential is performed on the early stage, and in- 
cluding the additional correlations is technically difficult. 
That is why we will derive the expression for the quan- 
tum correction in terms of classical probabilities, which 
are the solutions of the Liouville equation in a given po- 
tential. This result is important on its own, because it 
provides a tool for the description of the quantum effects 
in the ballistic cavities. The averaging, then, can be per- 
formed only on the final stage of the calculations. For the 
sake of concreteness, we consider two-dimensional case; 
generalization to the other dimensions is straightforward. 
We will omit the Planck constant in all the intermediate 
calculations. 



A. Introduction of basic quantities 

It is well known that transport coefficients can be cal- 
culated using the product of two exact Green functions 



Here G^^"*^ is the exact retarded (advanced) Green func- 
tion of the electron in the disordered potential U{r) and 
it satisfies the equation 

e ± zO - Hi] Gf '-^(ri, r2) = S{ri - rs), (3.2) 

where one-electron Hamiltonian is given by 

Hi = -^ + f/(ri). (3.3) 
For instance, the Kubo formula for the conductivity is 

(c^;ri,r2) 



de 

4m2 / 27r 



a/ 

de 



[^r" - [V^3 - V^J K,{u;; ri, r2, rg, r4) 



r4 = ri 
rs = i"2 



the expression for the polarization operator is 



n(t^;ri,r2) = s 



u5 (ri - r2) 
de df 



27r de 



K^{uj]ri,r2,r2,ri) 



(3.4) 



and so on. Here /(e) = (e(^^^)/^ + l) is the Fermi 
distribution function. Unfortunately, the exact calcula- 
tion of K is not possible and one has to resort on some 
approximations. 

In general, function iCe(a;; ri, r2, ra, r4) oscillates 
rapidly with the distance between its arguments. It con- 
tains non-oscillating part only if its arguments are paired: 
ri = r4, r2 = or, alternatively, ri = rs, r2 = r^. If 
they are not paired but still close to each other pairwise, 
then, it is very convenient to perform the Fourier trans- 
form over the difference of these close arguments: 



X.(..;ri,r2,r3,r4) = / ^^e-^(r.-r.)^.p.(r3-r.) 

(3.5) 



(27r)2 (27r)2 
xKf{u;; pi,Ri; P2,R2), 



Ri 



ri + Ti 



R2 



r2 + rs 



or, alternatively. 



i^.(..;ri,r2,r3,r4) = / .p.(..-.3),.p.(..-..) 
J [2nY (27r)^ 

xi^f (w; Pi,Ri; P2,R2), 

-I- r , 

(3.6) 



_ ri -I- ra _ rs + r4 
■tt-l — , -".2 — ■ 



Let us now derive the semiclassical equation for the 
function K'^ . From Eq. (|3^) and definition we can 
write the equation for function K in the form 
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Gf_. (rg, r4)<5(ri - r^) ~ Gf+. (n, r2)5(r3 ~ r4). (3.7) 

If the distance |ri — is much smaUer than the charac- 
teristic scale of the potential, we expand term H4 — Hi 
in Eq. ( ^3.71 ) in distance |ri — r^L and perform the Fourier 
transform analogous to Eq. ( |3.5D . The result can be ex- 
pressed in terms of the Liouvillean operator L: 



Li = 



dn d 



dn d 



9pi 9Ri 9Ri 9pi ' 
where 7i(p, r) is the Hamiltonian function 



W(p,r) = 



2m 



C/(r). 



(3.8) 



(3.9) 



With the help of Eqs. (^, (|3^) and (|3^), we obtain 



K^{uj- pi,Ri; p2,R2) 



(3.10) 



2^5 [e - 7^(P2, R2)] (27r)25(pi - p2)(5(Ri - R2 



Delta- functions in the right hand side of Eq. (3.10) should 



be understood in a sense of there subsequent convolution 
with a function smooth on a spatial scale larger than 
\f- When deriving Eq. (3.10), we used the semiclassical 
approximation for the Green functions 



Gf-4(ri,r2) 



dp 



oip(ri-r2 



(27r)2e-7^[p,(ri-f r2)/2]±z0 

(3.11) 



in the right-hand side of Eq. (|j) and neglected small 
frequency m in comparison with the large energy e ~ Ep- 
Liouvillean operator (^^) describes the motion of an 
electron in a stationary potential. Because the energy is 
conserved during such a motion, the function can be 
factorized to the form 

Kf {uj; pi,Ri; P2, R2) = X'e (t^; ni,Ri; n2,R2)x 



27r 



J[e-7^(pi,Ri)]5[e-7^(p2,R2)], 



(3.12) 



where diffuson I? is a smooth function of the electron en- 
ergy, n is the unit vector along the momentum direction, 
p = Pi?n = ny^2m [e — U{r)] and 1/ = m/27r is the den- 
sity of states. Diffuson I?e is the solution of the equation 



V^5i2, <5i2 = 2^(5(ni-n2)(5(Ri-R2). 

(3.13) 



It is important to emphasize that the diffuson 2? is a so- 
lution of the Liouville equation and not of the diffusion 
equation. In this sense, a more correct term for V is "Li- 
ouvillon" , however, we follow the terminology accepted 
in the theory of quantum disorder. 



Let us consider the classical chaotic motion such that 
the time of the randomization of momentum direction is 
finite. At small uj, which corresponds to the averaging 
over time scale much larger than the time of the momen- 
tum randomization, averaged over small region of its 
initial conditions, satisfies the diffusion equation 



V : 



1 



DV2 



(3.14) 



where D is the diffusion constant. The explicit relation of 
D to the characteristics of the potential U can be found 
in the limit of dilute scatterers Itr 3> a: in this limit the 
diffusion constant is given by D — VpTtr/2. It is worth 
emphasizing that Eq. (3.14) itself does not require such 
a small parameter, and it is always valid at large spatial 
scales and small frequencies . We will ignore the possible 
islands in the phase space isolated from the rest of the 
system. 

The semiclassical equation for function from 
Eq. (3.6) is found in a similar fashion: in the absence 
of the magnetic field and spin-orbit scattering it reads 



-iu + Li 



K^{uj; pi,Ri; P2,R2) = (3.15) 



2TrS [e - 7^(pi, Ri)] (27r)25(pi - p2)(5(Ri - R2). 
Function A'f can be factorized as 

K'^ (uj; pi,Ri; p2,R2)=C(cj; ni,Ri; n2,R2)x 



2tt 



(5[e-H(pi,Ri)](5[e-7^(p2,R2 



(3.16) 



Here Cooperon Ce is a smooth function of the electron 
energy satisfying the equation 



(3.17) 



Similar to the diffuson, the Cooperon, averaged over 
small region of its initial conditions, is a self-averaging 
quantity at large distances and small frequencies and 
in the absence of magnetic field and spin-orbit scatter- 
ing, it ca n be described by the expression analogous to 
Eq. (|t|): 



C = 



1 



DV2 



(3.18) 



B. Quantum corrections to classical probabilities. 

So far, we considered the lowest classical approxima- 
tion, in which the classical probabilities were determined 
by the deterministic equations of the first order. How- 
ever, the potential U contains not only the classical 
smooth part which is taken into account by the Liou- 
ville equation, but also the part, responsible for the small 
angle diffraction. Quantum weak localization correction 
originates from the interference of the diffracted electron 



6 



waves. The interference of the waves diffracted at differ- 
ent locations is added. It results, as we will show below, 
the quantum correction ceases to depend upon the details 
of the diffraction mechanism and becomes universal. The 
only quantity which depends on the diffraction angle is 
the time it takes to establish this universality. We will 
show, see also Sec. |l[ that the dependence of this time 
on the diffraction angle is only logarithmical. Therefore, 
with the logarithmic accuracy, we can include the effect 
of this diffraction into the classical Liouville equation by 
any convenient method, provided that we do it consis- 
tently for all the quantities and preserve the conservation 
of the number of particles. 

We will model the diffraction by adding the small 
amount of the quantum small angl e scatterers to the LHS 
of the Schrodinger equation (3.2). The effect of these 
scatterers will be twofold: 1) They will smoothen the 
sharp classical probabilities; 2) They will induce inter- 
action between the diffuson and Cooperon modes, that 
results in the weak localization correction. Finally, the 
strength and the density of these scatterers will be ad- 
justed so that the angle at which the classical probability 
is smeared during the travel to the distance a is equal to 



the genuine diffraction angle ^JXF /a. This procedure is 
legitimate because, as we already mentioned, the depen- 
dence of weak localization correction on the diffraction 
angle is only logarithmical. 

It is worth emphasizing, that even though weak local- 
ization correction takes its origin at the very short linear 
scale (ultraviolet cut-off) , the value of this correction at 
very large distances does not depend on this cut-off at 
all. Such phenomena are quite typical in physics, {e.g. 
in the theory of turbulence, theory of strong interaction, 
or in the Kondo effect). 

Let us now implement the procedure. Consider a single 
impurity located at the point s, and creating the poten- 
tial V{y) — Vb(s — r), so that the potential part of the 
Hamiltonian (3_^) is now given by U{y)^V{y). The char- 



acteristic size of this potential, d, is much larger than \p 
but much smaller than a. Our goal is to find the cor- 



rection to Eqs. (3.13) and (3.17) in the second order of 
perturbation theory in potential V . (Correction of the 
first order vanishes if T) and C are functions smooth on 
the spatial scale d.) In this order, correction to function 
(|3.lD has the form 



Ji^ (ri,r2,r3,r4) - j dr^re [G''ir,,v,)V{v5)G''{r,,r2)G^irs,Ve)V{ve)G^ive,u) + (3.19) 
G^(ri, r5)F(r5)G^(r5, r6)V^(r6)G«(r6, r2)G^(r3, r^) + G«(ri, r2)G^(r3, r5)V^(r5)G^(r5, r6)V^(r6)G«(r6, ra)] , 



where Green functions are the solutions of Eq. (3.2) with- 
out the impurity potential V. We will omit the energy 
arguments in the Green function, implying everywhere 
that the energies for the retarded and advanced Green 
functions are e + uj/2 and e — lu/2 respectively. 

In order to find the correction to the diffu son, wc con- 
sider the points ri,r4 and r2,r3 in Eq. ( 3.19 ) which 
are close to each other pa irwi se, perform the Fourier 
transfor m de fined by Eq. (3.5), and express the RHS 



We demonstrate the calculati on by evaluating the sec- 
ond term in the RHS of Eq. (3.19), let us denote it by 
SK2. 

Consider the product G-'^(ri, r5)G"^(r3, r4), points 
ri, r4 are close to each other, but points r5, are not. It 
means that for calculation of such a prod uct we can not 
use the semiclassical approximation ( 3.11 ) for the RHS of 
Eq. ( ^.7|) b ut sill can use the expansion (3.8) for the LHS 



of Eq. (3. IE) in terms of the diffusons and Cooperons. 



of Eq. (|3.7D. Solving Eq. with the help of Eq. ( ^3| ) 
we obtain 



G^(ri,r5)G^(r3,r4) = ^ J dv,v, j ^|EH_e'Pi(ri-.)-.p.(r.-r.)p 
[G«(r7, v,)5{v^ - r3) - G^{v:i,v^)5{v^ - v,)] 5 



ri + r4 r7 + rs 

2 ' ' 2 



(3.20) 



H P2, 



with H(p,r) being the Hamilton function (3.9). We 
will omit the frequency argument in the diffusons and 
Cooperons, implying everywhere that it equals to lj. 

We substitute Eq. ( ^.2C ) into the second term in the 



LHS of Eq. ( ^.1£ ) . We neglect the product of three re- 
tarded Green functions because this product is a strongly 
oscillating function of its arguments and vanishes after 
the averaging on a spatial scale larger than \p. The 
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remaining product G'^(rg, r2)G^ (r^^s) is approximated product of two retarded Green functions and performing 
by the expression similar to Eq. ( p. 20 ) because points r2 the Fourier transform over the differences ri — r4 and 
and rs are close to each other. Neglecting, once again, the r2 — rs, we find 



<^2?2(1,2) 



SK2 (r+,r+,r2 



(ir3(ir4(iR3(iR4 



27riy 



dps dp4 
(27r)2 (27r)2 
drii dn2 



g_,p3r3+»p4r4p 3) p (4. 2) qR (r+, r+)G^(r4 , r. 



(3.21) 



(27r) (2^) 

Here we introduced the short hand notation j = (nj,Rj) and = Rj ± 

What remains is to find the semiclassical expression for the product G^G^ in Eq. ( 3.2l| ). We no tice that the points 
rj]',r^ lie within the radius of the potential V{r). In order for the product G^G^ in Eq. ( |3.21 ) not to vanish, the 
points , must be close to the points r;^,r4 . Because all the four points are close to each other, one can write, 
cf. Eq. (p), 

G^(r+,r+)G'^(r4 ,r3-) = J dn^dns 9 [n4(r+ - r+)] 6 [nz{rl ~ v^)] e'P^^^{4-^t)+^P^^^i^^ -^',) + 



V 

2^ 



n I "^5, 



(3.22) 



Here, the first term is the explicitly separated contri- 
bution of the short straight line trajectories connecting 
points ri,r2 and r3,r4. These short trajectories can be 
well described by the Cooperon or by the diffuson. The 
second term describes the contribution of all the other 
trajectories connecting these points. It can be shown by 
explicit calculation that the representation of Eq. ( 3.22| ) 
in terms of the diffuson only would lead to the loss of this 
second term. This is because the Cooperon describes the 
interference effects corresponding to the oscillating part 
of the di ffuso n which is lost in the semiclassical approxi- 
mation ( 3.12 ). 

Now, we are ready to find the correction coming from 
the singl e qua ntum scatterer. We substitute Eq. (3.22) 
into Eq. ( 3.21 ) and perform the integration while neglect- 
ing the dependence of the diffusons and Cooperon on 
their spatial coordinates on the scale of the order of the 
scatterer size. We consider the remaining two terms in 



Eq. ( 3.19 ) in a similar manner. The overall result is 



8V = 5V 



St 



6V' 



(3.23) 



,5X'^*(1,2) = J d3d4Ps(3,4)r'(l,3)[2?(4,2)-X>(3,2)], 



zni' 

[I?(1,3)-I?(1,4)][I?(3,2) 



I?(4,2)] 



Here we use the short hand notation j = (rij , Rj ) , inte- 
gration over the phase space on the energy shell is de- 
fined as dj = dnjdHj /2Tr, the time reversed coordinate j 
is given by j = (— n^, Rj), and the kernel V describing 
the scattering by an impurity is 



( 1 , 2) = 27ri/(5 (s - Ri ) (5(s - R2 



The first term SD^* in Eq. ( 3.23 ) coincides with that 
obtained for otherwise free moving electrons. The second 
term, SD^ , describes the interference effect arising be- 
cause the chaotically moving in classical potential U{r) 
electron may return to the vicinity of the impurity one 
more time. 

The correction to the Cooperon due to the single im- 
purity can be obtained from Eq. (3.19) by considering 



close pairs ri, and ^2, r4; this results in the expression 
similar to Eq. (3.23) with the replacement T) ^ C. 

So far, we considered the correction due to a sin- 
gle weak impurity. If the number of these impurities 
is large, we can, in the lowest approximation, consider 
the contributions from the different impurities indepen- 
dently of each other, by the substitution to the RHS of 
Eqs. (3.23) of the diffusons and Cooperons renormalized 
by all the other impurities. As the result, we arrive to the 
Boltzmann-like equations for the diffuson and Cooperon 



—iuj + Li 
—iu) + Li 



V{1, 2) = Si2+J2j ^3^,(1, 3) |[P(3, 2) - V{1, 2)] + [I?(i, 2) - I?(3, 2)] + 
C(l,2)=^i2 + ^ /'d37'.(l,3)|[C(3,2)-C(l,2)]-f [C(T,2)-C(3,2)] 



27rzy 

I?(l,3) -f I?(3,l) 



(3.24a) 
(3.24b) 



where the nota tion for the coordinates j, j was intro- 
duce d aft er Eq. ( 3.23| ), and the S - symbol was defined in 
Eq. ( ^I3|) . 

Assuming that the distribution of the quantum impu- 
rities is uniform with the density n^, we can make the 
continuous appro xima tion and replace —^rii J ds in 
the RHS of Eqs. ( 3.24 ). Finally, taking into ac coun t that 
the scattering angle is small , we r educe Eqs. ( |3.24| ) to a 
differential form. Equation (3.24a) becomes 



-lUJ - 



Li 



1 ^ 

Tg d4>l 



V{1) 



d c(i,i) d 
5i2 - m-^O-) 

0<Pl HVTq 0(pi 



(3.25a) 

Here, angle is defined so that rij — (cos0j, sin cjjj), the 
notation for the coordi nates j = (nj,Rj), j = (— nj,Rj), 
is the same as in Eq. ( ^.23 ), and the S - symbol was de- 
fined in Eq. (3.13). The second argument is the same 
for all the diffusons in Eq. ( 3.25aD and that is why we 
omitted it. Analogously, Eq. (3.24b) reduces to 
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(3.25b) 



The se cond a rgument is the same for all the Cooperons 
in Eq. ( 3.25b ) and it is omitted. Quantum transport life 
time in Eqs. ( 3.25| ) is given by 



— = 2'Kvni 



271 2 



Equations (3.25) describe how the classical Liouville 
equation changes under the effect of the small angle scat- 
tering (diffraction). We see that the quantum effects re- 
sult in two contributions to the Liouville equation. First 
contribution provides the angular diffusion and, thus, it 
leads to the smearing of the sharp classical probabilities. 
Usually, for the calculation of the transport coefficients, 
such as the diffusion constant or the conductivity, the av- 
eraging over the initial and final coordinates is performed 
anyway. Therefore, the angular diffusion itself provides 
only negligible correction to the classical transport coef- 
ficients which are controlled by classical potential U . On 
the contrary, the second contribution givin g the q uantum 
correction [last terms in the RHS of Eq. ( 3.25a )] is pro- 
portional to the classical probability C(T,1) where the 
initial and finite points of the phase space are related 
by the time inversion. In the absence of the spreading 
due to the angular diffusion, Tg ^ oo, this probability 
vanishes identically, see Sec. |l[ In order to obtain the 
correction at finite time (or finite frequency), one must 
keep Tq finite even in the final results. 

Let us estimate the value one should ascribe to Tq 
for the description of the diffraction effects in the sys- 
tem. As we already discussed, for the calculation with 



the logarithmic accuracy, we do not need the numeri- 
cal coefficient. The parametric dependence of Tq can be 
established by using the following argument. Consider 
two independent electrons, starting with the same ini- 
tial conditions. If there were no diffraction, they would 
propagate together forever. Due to the angular diffusion 
(diffraction), the directions of these trajectories deviates 
first oc y/t and then exponentially, ^^^-^ « 2\{5(j)'^) + i^ , 
where angle 5(j) stands for the angle between the momenta 
of two electrons, and A is the Lyapunov exponent. It 
yields {5(j)'^it)) « (4Ar,)-i (e^^* - l). Thus, the charac- 
teristic time during which the angular diffusion switches 
to the exponential growth is always — On 
the other hand, quantum spreading of the wave packet 
during this time interval is given by 5x^ ~ Xpypte- 
Taking into account the relation 6x ~ (jft^vp, we find 
te/Tq ~ Xp/{vpte). It yields the esti mate f or the quan- 



tum transport time entering into Eqs. (3.25) correspond- 
ing to the small angle diffraction 



1 

Tn 



A 



2 Af 
Vp ' 



(3.26) 



It is important to emphasize that the very same Tq en- 
ters into the angular diffusion term and into the diffuson- 
Cooperon interaction. This circumstance is extremely 
crucial for the universality of the quantum correction at 
large time (w ^ 0), even though parameter Tq itself does 
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not enter into the result, see Sec. 

Let us now turn to the calculation of the lowest quan- 
tum correction to the diffuson. Taking into account the 
last term in the RHS of Eq. (3.25a) in the first order of 
perturbation theory, we obtain 

X>(1,2) =X>°(1,2) -f AX>(1,2); (3.27a) 
C0(3,3) ai?0(l,3) 9X»0(3,2) 



^5(1, 2) ^ J d3 



-lu; + Li ^ 



— ILU + Ll — 



m/Ti 

Tq d(l)l 
1 92 1 



d(t>'i 



d(t>'i 



Tq d(t)l 



2?°(l,2)-(5i2; 
C°(l,2)==<5i2, 



(3.27b) 
(3.27c) 
(3.27d) 



where j = (nj,Rj), integration over the phase space on 
the energy shell is defined as dj = dnjd'R.j /2tt, the time 
reversed coordinate j is given by 7 = (— nj,Rj), and the 
d - symbol was defined in Eq. (PT3|) . 

Equation ( |3.27b| ) can be rewritten in a different form. 
Even though more lengthy than Eq. (3.271), this form 
turns out to be more convenient for further applications: 

AI?(1,2) =I?°(l,2)^^^^ + ^^^Mp°(l,2) (3.28) 
2t:i' 2'kv 



d3X'°(l,3)X'°(3,2) 
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2iu} — iyj H 



C°(3,3) 

2-KV 
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In order to derive Eq. ( |3.28D from Eq. ( |3.27bD , we sub- 
tracted from the RHS of Eq. ( 3.28 ) the expression 



C(3,3) 
27rv 



which vanishes because integrand is the total deriva- 
tive along the classical trajectory. Then, we int egrate d 
Eq. ( ^.27b| ) by pa rts an d, with the help of Eq. ( |3.27cD , 
we arrived to Eq. (^.28[ ). 

Equations ( |3.27b| ) and ( |3.28| ) are the main results of 
this section. They give the value of the lowest quantum 
correction to the classical correlator in terms of the non- 
averaged solutions of the Liouville equation (with small 
angular diffraction added) for a given system. Besides the 
found correction, there exist the other corrections [e.g. 
from the hi gher terms in the expansion ( |3.8| )], however, 
Eqs. (3.27b) and (3.28) are dominant at low frequencies. 
The quantu m m esosco pic fluctuations are neglected in 
Eqs. ( 3.27b ) and ( 3.28 ), which implies either the temper- 
ature is high enough or the averaging over the position of 
the Fermi level is performed. Then, if the relevant time 
and spatial scales are large, the quantum correction be- 
comes a self-averaging quantity expression for which will 
be obtained in the next section. 



IV. AVERAGED QUANTUM CORRECTIONS 

We will consider the quantum correction at large dis- 
tance and time scales. In this case, the classical proba- 
bility does not depend on the direction of the momentum 
and it is given by Eq. ( 3.14| ). Our goal now is to find the 
expression for the quantum correction in the same ap- 
proximation. We will bear in mind the systems in which 
the diffusion constant is large enough, D/avp > 1. It is 
the case for the antidot arrays. The conductance of the 
net of the ballistic cavities requires a separate considera- 
tion. 

For the calculation we use Eq. ( 3.28| ). While per- 
forming the averaging, we make use of the fact that the 
Cooperon part of the expression can be averaged inde- 
pendently on the diffuson part. This is because the clas- 
sical trajectories corresponding to these quantities lie es- 
sentially in the different spatial regions, (see e.g. Fig. 
where segments iB and fB correspond to diffusons and 
segment BEE corresponds to the Cooperon) and, there- 
fore, they are governed by the different potentials and are 
not correl ated. Performing such an averaging, we obtain 
from Eq. ( |3^ ): 



AX>(1,2) 



(pO(l,2)) + (I?"(l,2))- 



2iuj / d3(X'°(l,3)X'°(3,2)) 



(4.1) 



2lTV 



where (...) stands for the averaging either over the real- 
ization of potential U or over the position of the "center 



of mass" of the Cooperon and diffuson. The last two 
terms in brackets in Eq. ( 3.28 ) vanish after averaging 
because the averaged cooperon does not depend on the 
coordinates ns, R3. 

On the other hand, as we have already explained in 
Sec. H, the correlations in the motion of both ends of 
the Cooperon can not be neglected. The same is also 
true about the correlation betwe en m otion of the ends 3 
and 3 in the third term of Eq. ( ^.28| ). In what follows, 
we will separate the description of the problem into the 
Lyapun ov and diff usion regions. It will be done in sub- 
sections IV A and [VB for the Cooperon and diffusons 
respectively, and the resulting cor rectio n to the conduc- 
tivity will be found in subsection [V C . The description 



of t he Ly apunov region is presented in subsections IVD 
and lIVEl. 



A. Cooperon in the diffusive and Lyapunov regions 

In order to find (C°(l,l)) we consider more general 
quantity C(0, p) defined as 

C{^,P) = U ^C"(n+,R-;-n-,R+), (4.2) 



2tt 



n cos I ± [n X 1^] sin |, R± = R ± ^ [n± x 1^] , 

where S is the area of the sample, and Iz is the unit 
vector perpendicular to the plane. Function C(0, 0) co- 
incides with the necessary quantity (C°(l, 1)). 

It is easy to find C((/), p) in the diffusion region. At 
Pfthr, it is given by 



C(0,p) = 



1 



-iuj- D^l' 



(4.3) 



At p < y^D/oj the Cooperon depends only logarithmi- 
cally on p and at a ^ p ^ ^tr, it becomes independent of 
p. With the logarithmic accuracy, we have 



C(0,p) 



In I ) , p > a. 



(4.4) 



Equation (4.4) serves as the boundary condition for 
C (0, p) at the boundary between the diffusive and Lya- 
punov regions: 



b,p = a signc; 



1 



A-kD 



■In 



1 



UJTtr 



(4.5) 



Meaning of Eq. (4.5) is that both ends of the Cooperon 
enter into the Lyapunov region with the random mo- 
menta, and thus the probability of this entrance is given 
by the solution of the diffusion equation. 

The next step is to find C (0, p) in the Lyapunov re- 
gion. We add to Eq. (3.27d), the equation conjugate to 
it, which gives 



10 
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1 (9 



2 n 



C"(l,2) = 2(5i2. 

(4.6) 



Formula ( |4.q ) enables us to find the equation for qua ntit y 
C° {n,R;<j),p) = (n+,R-;-n-,R+) from Eq. Q. 
Expanding potential U up to the first order in p, and 
using the fact that the angle (j) is small, we obtain 



-2iuj + Lc - 
Here operator 
L = 



2t„ dS 



C"(n,R;0,p) = O. (4.7) 



9R 



dUjH) d_ 
(9R ■ 9P 



(4.8) 



describes the motion of the "center of mass" of the 
Cooperon along a classical trajectory and operator Lj. 
characterizes how the distance between the ends changes 
in a course of this motion: 



-Vf4> 



d 



dp pFdR\ 



d_ 



(4.9) 



with R± being the projection of R onto the direction per- 
pendicular to n. In Eq. we neglected the effect of 
the angular diffusion on the motion of the center of mass 
because the averaging over the posi tion of the center of 
mass n, R is performed in Eq. (4.2) anyway. 

Now, we have to find function C{p, (f) in the Lya- 
puno v re gion, satisfying the bounda ry c onditi on g iven by 
Eq. (HH) and consistent with Eqs. (^Jj) and (^). Solu- 
tion can be represented in a compact form analogous to 
Eq. (U) 



C(0,p) = !filfi^ln 



1 



Function w [lo] (j), p) is defined as 



u; (w; 0, p) = - / VF(w;n,R;0,p), 



(4.10) 



(4.11) 



where S is the area of the sample and W is the solution 
of the equation 



-2iijj + Lc + Lr ~ - — 
2t„ 
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M/(c^;n,R;0,p) = 0, 

(4.12) 

supplied with the boundary condition 

W^(tj;n,R;0,p = a sign0) = 1. (4.13) 

The necessary quantit y (C° (l,i)) is, thus, found by 
putting p,(j) = in Eq. ([l.lOD 



(C"(l,l)) = !^il^ln^ 1 



AttD 



\UJTtr 



(4.14) 



B. Diffusons in the diffusive and Lyapunov regions. 



In this subsection wc find the average 



/ d3(X>°(l,3)2:'°(3,2)) entering into Eq. (Q. We use 
the procedure similar to the calculation of the Cooperon 
in the previous subsection. We consider more general 
quantities M, A4 defined as 

7W(l,2;n,R; 0, p) = I?° (l; -n", R+ ) I?°(n+, R"; 2) 
M(l,2;^,p)=y'^(M(l,2;n,R;<^,p)), (4.15) 

where the coordinates n*,R^ are defined in Eq. (|4.2|). 
Function M (1, 2; 0, 0) coincides with the necessary quan- 
tity /fi3(2?0(i,3)2?0(3 3)). 

In the diffusive region p^ a two diffusons are governed 
by the different potentials and, therefore, can be av eraged 
independently; each of them is given by Eq. (3.14). Fur- 
thermore, if p ^ \J~DJuj^ function Af(l,2;0, p) becomes 
independent of p, and it is given by 



M(l,2;0,p)= / d3(2?"(l,3))(2?°(3,2)). 



(4.16) 



Equation ( 4.16 ) serves as the boundary condition for 
Af (1, 2; (/), p) at the boundary between the diffusion and 
Lyapunov regions: 

M(l,2;(/.,p = asign0)= /d3(I?°(l, 3))(I?"(3, 2)). (4.17) 



Meaning of Eq. ( 4.17 ) is that the ends of both diffu- 
sons enter into the Lyapunov region with the random 
momenta. 

The next step is to find M(l, 2; (/), p) in the Lyapunov 
region. It follows from Eq. ( |3.27c ), that the product of 
two diffusons I?'^(l; 3) I?°(4; 2) satisfies the equation 



-2ia; + L-i 
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%I?°(4;2)-f 5242?"(1;3) 



X'°(l;3) X>"(4;2) 



(4.18) 



Equation ( |4.18| ) e nables us to find the equation for quan- 
tity M. from Eq. (4.15). We expand the potential U up 
to the first order in p, and use the fact that the angle 
is small. This yields 



-2iu} + Lc + Lr 

2t„ 
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X(l,2;n,R; 0,p) 



2n5 (m + n-) (5 (Ri - R+) P"(n+, R"; 2) + 

27r(5 (n2 - n+) 5 (R2 - R") -n", R+) , (4.19) 



where the operators Lc and Lj. are defined in Eqs. (4.S) 
and (4.9) respectively. In Eq. (4.19), we neglected the 
effect of the angular diffusion on the motion of the center 
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of mass because the averaging over the posit ion of the 
center of mass n, R is performed in Eq. ( 4.15| ). 

We have to find function M{1,2; p, (p) in the Lya- 
punov region , satisfying the boundary con dition give n 
by Eq. ( 4.17 ) and consistent with Eqs. ( 4.15 ) and ( 4.19 ). 
We represent functions M, Ai as the sum of two terms 
M = Mi+ Ma, M = M1+M2, 



M,(l, 2; 0, p) = J^^{M.{1, 2; n, R; </>, p)), (4.20) 

for j = 1, 2. Function AI1 is a solution of the inhomoge- 
neous equation 



-2iuj + Lc + Lr 



2t„ c'(/)2 



Xi(l,2;n,R;0,p) 



2ttS (m +n-)(5(Ri -R+)l?°(n+,R-;2) + 

2ttS (na - n+) S (R2 - R") P°(l; -n", R+) , (4.21) 

without any boundary conditions imposed and function 
is the solution of the homogeneous equation 



-2iuj + Lc + Lr 
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2t, 902 



A^2(l,2;n,R; </),p) = 0, 
(4.22) 



with the boundary condition 



M2 (1,2; <j),p = a sign0) = yd3(I?°(l, 3))(P°(3, 2)) - 

Mi(l,2;(/),p = asign0). (4.23) 

First , we find function Mi. We integrate both sides of 
Eq. ( 4.21 ) over R, n and average them. This gives 



-2iLU — 



2t„ dd) 



A/i(l,2;p, 



(4.24) 



dndR 
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{LrMi (1, 2; n, R; 0, p)) = (I?"(i; 2)) + (I?°(l; 2)). 



Calculating the RHS of Eq. ( |4.24| ), we neglect p < a < 

V Duj in the arguments of the averaged diffusons. Right 
hand side of Eq. (4.24) is independent on p and (f>. There- 
fore, we can seek for the function Mi(p, </>) al so in depen- 
dent of p, (j). The last term in the LHS of Eq. ( 4.24 ), then, 
vanishes and we obtain 



Mi(l,2;p, 



(2?°(1;2)) + (I?°(1;2)) 
—2iuj 



(4.25) 



Substituting Eq. ( [4.25| ) into Eq. ( [4.23D , we find the 
boundary condition for the function M2 



M2(l,2;0,p = a sign0) = j d3(I?°(l, 3)) (I?"(3, 2)) 
(I?0(T;2)) + (I?"(1;2)) 



-210; 



(4.26) 



Equa tion (4.22), supplie d w ith th e bo undary condition 
( 4.26 ), is similar to Eqs. (4/7) and ( |4.5| ) for the Cooperon 
consi dered in the previous subsection. Thus, we use 
Eq. ([4.IOD to obtain 



M2 [1,2] (I), p) = w{uj](t),p) 



d3{V°{l,3)){V°{3,2))~ 



(I?"(I;2)) + (I?°(1;2)) 
—2iu! 



(4.27) 



where function w is defined by Eq. (4.11). 

The necessary quantity / (i3(r>°(l,3)X>"(3 2)) is, thus, 
found by summing the contributions (4.25) and (4.27) 
and putting p,(f) = 0. We obtain 



'd3(P°(l, 3)I?"(3, 2))^w{uj; 0, 0) /d3(I?°(l, 3))(I?"(3, 2)) 
l-w(w;0,0) 



-2iuJ 



[(I?«(i;2)) + (I?0(1;2))] . (4.28) 



C. Quantum correction to the conductivity. 



Now, we are prepared to find the correct ion t o the 
conductivity. Substituting Eqs. (4.14) and (4.28) into 
Eq. (O), and using Eq. (|3.14D for (P"), we find 



u;2(^;0,0)lnM- 



(4.29) 



wher e fun ction w is give n by Eq. ( 4.11 ). Comparing 
Eq. ( 4.29| ) with Eq. (3.14), we see that all the quantum 
correction can be ascribed to the change AD in the diffu- 
sion constant. Restoring the Planck constant, we obtain 



AD{lu) 



(4.30) 



The correction to the conductivity Act is related to the 
correction AV by Einstein relation Act = se'^vAD, where 
s = 2 is the spin degeneracy. We immediately find 



Act 



u;2(a;;0,0)ln 



/ 1 



(4.31) 



Comparing Eq. (4.31) with Eq. (1.1), we obtain the renor- 
malization function T{lo): 



r(w) = ^2 (w;0,0). 



(4.32) 



Here function w is defined by Eq. (4.11) 
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D. Universality of the weak localization correction 
at u — > 0. 



The universality of weak localization correction at low 
frequencies, r(0) ~ 1 can be proven imm ediate ly. In- 
deed, function = 1 is a solution of Eq. ( 4.12 ) and it 
satisfies the boundary condition W{p = a) = 1. Because 
W = 1 is the solution of nonaveraged equation in specific 
disordered potential, the averaged func tion w also equal s 
to unity. Then, it follows from Eqs. ( fl.ll| ) and ( [4.32[ ), 
that r(0) = 1, which completes the proof of the univer- 
sality. This fact is well-known for the weak short range 
disorder, where the Born approximation applies. We are 
not aware of any proof of the universality for the disorder 
of the arbitrary strength and the spatial scale. 

We emphasize that the proof did not imply any small 
classical parameters in the problem, and it requires 
only the applicability of the semiclassical approximation, 
Xp <^ a,ltr- Universality is based on two elements: 1) 
the conservation of the total number of particles on all 
the spatial and time scales and 2) existence of a diffusive 
motion at large spatial and time scales. Both these facts 
depend neither on the strength of the scatterers nor on 
their spatial size. 

It is worth ment ioni ng also that the upper cut-off of the 
logarithm in Eq. ( |l.l| ) is determined by purely classical 
quantity Ttr and does not contain Ehrenfest time as one 
could expect. This result is due to the fact that the both 
lower and upper limit of the logarithm in the solution of 
the diffusion equation are related to the spatial scale and 
not to the time scale. The upper limit of the logarithm 



y^D/uj is the typical distance at which the electron can 
diffuse during time ~ l/w. The lower linear scale is the 
largest of two distances; 1) the distance between the ini- 
tial and final points, or 2) the transport mean free path 
- smallest scale at which the diffusion approximation is 
applicable. Because, for the problem in the diffusive re- 
gion, we are interested in the probability for an electron 
to approach its starting point at the distance of the or- 
der of a < Itr, (and by no means y/DtE), we have to 
use Itr as the short distance cutoff. It immediately gives 

ln(yW^/^tr) = ln(l/V^)- 

Thus, we conclude that the weak localization correc- 
tion has precisely the same universal form as in the quan- 
tum chaos regime. However, unlike in the QD regime, 
this universality persists only up to some frequency which 
is much smaller than Ttr and breaks down at larger fre- 
quencies. The description of such a breakdown is a sub- 
ject of the following subsection. 



E. Ehrenfest time and T{lu) at finite frequency. 



Our goal now is to find w at frequencies t^l < w < 
Tj~^. We would like to show that the functional form of 
w is log-normal even if the parameter a/kr is not small, 



and derivation of the equation analogous to the Boltz- 
mann kinetic equation is not possi ble. L et us, first, ne- 
glect the angular diffusion in Eq. (4.12) at all, we will 
take it into account in the end of the subsection. We 



rewrite Eq. (4.12) in the time representation 



d_ 

dt 



vf4> 



d 



dp ppdH 



d__mj_ d_ 

d^U d 
— tPjt. 



iy(i;n,R; 0,p) = 0, 



duj 
2^' 



*VF(w) 



(4.33) 



where we used the explicit form of operators Lc,r from 
Eqs. (4.8) and (4.9). Then, we separate the motion of 
the center of mass and the relative motion of the ends of 
the Cooperon. Namely, we factorize function W as 

— {t; no, Ro; (f>, p) x 

5 [R - R(t, Ro, no)] 5 [n - n(t, Ro, no)] , (4.34) 

where the trajectory of the center of mass R(t),n(i) is 
found from the classical equations of motion 



aR' 



R 



P 

m ' 



n{t) 



|P(i)l' 

R(0)=Ro; P(0) = noPF(Ro) 
and function W±_ obeys the equation 



«F(i) =«F[R(i,Ro)], F{t) 



VF_L =0, 



(4.35) 



(4.36) 



ppdRi 



R, — Fi,(t ,R,o) 



Equation ( 4.36| ) is invariant with respect to the scale 
transformation of variables p and 0. It invites to intro- 
duce the new variables 



In 



(- 

\a 



arctan — . (4.37) 
P 



Upon this substitution, Eq. (4.36) takes the form 



|-5,(.)sin(2.)A_ 



[Si(t)cos(2a) + B2{t)] — \W^^0- 



Za 2 



(4.38) 



Formal solution of Eq. (4.38) is (we omit arguments 
no ; Ro hereinafter) 



W± {t; z, a) — exp 



W± [0-z,aoia,t)] 



B3{t,a)= [ dtiBi{ti)sm{2a[ao{a,t) ,ti]}, (4.39) 
Jo 
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where function a{ao,t) satisfies the equation of motion 
da 



dt 



^ Bi{t)cos{2a) + B2{t), a(ao,0)=ao, (4.40) 



and function cyo{t, a) is imphcitly defined by the relation 
a [ao{t, a), t] = a. (4-41) 



Equation ( 4.39| ) enables us to find the time evolution of 
func tion w{t) from Eq. (4.11). Indeed, substitution of 
Eq. ( 4.34 ) into Eq. ( 4.11 ) immediately yields 



w{t;(t>,p) 



27r5 



VFi(i;no,Ro;0,p). (4.42) 



The time dependence of the function W± is given by 
Eq. (4.39); using this formula we obtain 



w (t; z, a) : 



dRodno 



-exp 



w[0; z, ao(i, a)] . 

(4.43) 



We are interested in the time dynamics of the system 
at time t much larger than Ttr- At such large times, func- 
tion a{aQ,t) averaged over an arbitrary small region of 
R(),nQ is a self-averaging quantity and it no longer de- 
pends on the initial condition uq. (This fact is similar 
to the randomization of the direction of momentum in 
the derivation of the diffusion equation). Therefore, the 
function i?3 from Eq. (4.39) becomes independent of a. 
Thus, at large times w (t\ z, a) is also independent of a 
and its evolution is governed by the Focker-Planck type 
equation: 



d_ 

dt 



w{t, z) — 0, 



(4.44) 



where J-{x) is defined as 



Bit) 



lini — In 

t^oo t 



dRodno 

2ttS 



exp [xB{t)] 



dtBi(t) sin [2a(ao,t)] 



(4.45) 



In Eq. (4.45), the initial condition ao may be chosen ar- 
bitrary. Furthermore, we will need function w at large 
times. In this case is a smooth function on z, and we 
expand T in the Taylor series: 



T{x) = Xx 



X = lim — 

t^oc t 



lim — 

t^oo t 



As 



X2X' 

2 

dRodno 

dRodno 
2ttS 



(4.46) 



xh^ 



Returning to the frequency representation, we obtain 
the equation describing the drift and diffusion of the log- 
arithms of the coordinates: 



.9 X2 

-2iuj — X— 



82 



2 dz 



wiuj] z) ~ 0. 



(4.47) 



With the same accuracy, the boundary conditions 



Eq. ( [4.13D take the form 

w{uj] z = 0) = 1 



(4.48) 



For a generic system the actual calculation of the coef- 
ficients A, A2 can be performed, e. g. by the numerical 
study of the system of equations ( 4.35| ) and (4.40) at 
times of the order of Ttr and then using Eq. (|4.46|). An- 
alytic calculation of coefficients A, A2 requires additional 
model assumptions. Outline of such calculation for the 
weak smooth disorder i s pre sented in Appendix |^. 

The solution of E q. (4.47 ) at uoTtr ^ 1 and with the 
boundary condition ( |4.4g ) has the form 



w{uj; z) 



exp 



A 



A3 



(4.49) 



However, in order to find the renormalization function 

r(ti;), we need to know w{p^(f) — 0), see Eq. ( 4.32 )^ It 

corresponds to taking the limit z —00 in Eq. (4.49). 
One immediately realizes that w{p, cf) = Q) ~ Q at any 
finite frequency w, which would mean that the time it 
takes for the quantum correction to reach its universal 
value is infinite. The reason for this unphysical r esult lies 
in neglecting the angular diffusion term in Eq. ( 4.33 ). It 
is this term that is responsible the quantum spreading of 
the classical probability and it makes the Ehrenfest time 
finite. 



In terms of the variables ( 4.37 ), the angular diffusion 
operator is given by 



52 



-2z 



92 

z 



cos2q; — e 
oz 



-2z 



d_ 

dz 



O 



d_ 

da 

(4.50) 



Because function w is independent of a, we can neglect 
all the terms O at all. Furthermore, the condition 
Atq ^ 1 enables us to consider the angular diffusion 



(E 



50) in the l owest order of perturbation theory. As 
the result, Eq. (1.47) acquires the form 



2iL0 



d_ A2 d'^ 
'd'z ^ Yd^ 



d /I — 7 9 



2Tg dz 



2 dz 



+7 



w = 0, 
(4.51) 



where the numerical coefficient 7 ^ 1 is given by 



7 = lim — dt 

t^OO t Jq 



dRodno 
2ttS 



cos [2Q;(ao, t)] 



We now solve Eq. (4.51) with the logarithmical accu- 
racy, taking into accoimt the condition Ar^ ^ 1. The 
result is 
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At |zj < \nXTg, expression ( |4.52D matches with The factor e"*^/'^'^ in Eq. (O) can be easily understood 



Eq. (4.49) 



By taking the limit z —> — oo in Eq. (4.52) and making 
use of Eq. (4.32), we obtai n Eq (1.2) with = j In At,. 
Finally, we use estimate ( 3.26| ), replace with the loga- 
rithmic accuracy vp/X to the cha racteristic size of the 
potential a, and arrive to Eq. (1.3). 



V. RELEVANT PERTURBATIONS. 

So far, we considered only the frequency dependence 
of weak localization correction in the quantum chaos. In 
this section we concentrate on two more factors which 
affect our results: 1) finite phase relaxation time t^; 2) 
presence of the magnetic field; 

A. Effect of finite phase relaxation time r^p 

As it was discussed in Sccjl|, the weak localization 
correction has its origin in the interference between the 
coherent classical paths. If the particle experiences 
the inelastic scattering during its motion, this coher- 
ence is destjjOjjcd and the weak localization correction is 
suppressedQli30. This effect is described conventionally 
by the introduction of the phase relaxation time t^^,, (see 
Ref. [ |l^ for a lucid discussion of the physical meanin g 



of T^), into the Liouville equation for Cooperon (3.27d) 



1 ^ 1 92 
— +Li — . 

Tq 0<p\ 



c = s 



12- 



(5.1) 



The equation for the diff uson ( 3.27c) remains unchanged 
as well as the relations ( 0.27b ) and ( 3.28 ) between the 
correction to the classical probability and the Cooperon 
and diffusons. 

Thus, we have to modify the Cooperon part of 



Eq. (4.31); namely Eq. (4.10) acquires the form 



C(0,p) = ^(" + ^/"^''^'^^ In 



(5.2) 



Comparing Eqs . ( 4.14) and ( p. 21) , we obtain with the 
help of Eqs. and ([4.32D 
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r (w) r w -t- — 



1/2 



In 



A relevant trajectory may close not earlier than it leaves 



the Lyapunov region; factor e 



-tElt 



is nothing but the 



probability for an electron not to be scattered inelasti- 
cally while it is in the Lyapunov region. Let us notice also 
that the dependence of the weak localization correction 
on the phase relaxation time is always slower than an ex- 
ponential. The reason for this is the following. The prob- 
ability for a trajectory to leave the Lyapunov region dur- 
ing time interval Ti^/2 is determined by the correspond- 
ing Lyapunov exponent and, thus, it can be increased 
due to the fluctuation of this exponent. The probability 
to find such a fluctuation is given by the Gaussian dis- 
tribution. The optimization of the product of these two 
probabilities immediately yields the exponential factor in 



Eq. ( pD . 

At this point, we should caution the reader, that the 
fact that the same T^p enters into the logarithmic fac- 



tor and into the renormalization factor F in Eq. ( |5.3| ) 
is somewhat model dependent. Strictly speaking, this 
statement is valid only if the phase breaking occurs via 
single inelastic process with the large energy transfer. If 
the main mechanism of the phase breaking is associated 
with the large mimber of scattering events with the small 
energy transferli3'tZl, the phase breaking occurs when the 
distance p between the Cooperon ends is large enough, 
\/~DJT ^ p ^ ^Dt^. Thus, this mechanism does not af- 
fect the Cooperon in the Lyapunov region at all. Further 
discussion of the microscopic mechanisms of the phase 
breaking is beyond the scope of the present paper. 



B. Effect of magnetic field 

Similar to the phase relaxation time, the effect of the 
magnetic field on the weak localization correction is taken 
into account by thci-daaagc in the equation of motion for 
the Cooperon only:Qli3't2l 



—ill) + L\ 



2ie^ A ^ ' 

c ^ ^ Tqd(l)l 



C^6 



12, 



(5.5) 



where Ai = A(Ri) is the vector potential of the exter- 
nal magnetic field. Cooperon given by Eq. (5.5) is not 
a gauge invariant quantity but C(l, 1) is. It is very con- 
venient to separate the gauge noninvariant part of the 
Cooperon explicitly by writing 



C = exp ^— J Adr^ Cgi, 



(5.6) 



(5.3) 



For Lu = and 3> Tf,., expression (5.3) acquires the 
form 



where integration in the first factor is carried out along 
the straight line connecting the Cooperon ends. Sub- 
stituting Eq. ( ^ ) into Eq. (5.5), we obtain the gauge 
invariant part of the Cooperon 
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Li 



[z X r] 



12, 



(5.7) 



1/2 

where r = Ri — Ri and \h = {ch/eH) ' is the mag- 
netic length. When the ends of the Cooperon coincide, 
Cgi — C, and, therefore, the correction to the conductivity 



(4.31) is modified as 



Aa = 



se 



w{uj;0,0){{C^i{l,l))D). 



(5.8) 



Our purpose now is to obtain the expression for Cgi. 
Similar to the case of zero magnetic field, we would like to 
separate the problem into Lyapunov and diffusion region. 
This separation, however, is valid only if the condition 



Ah > k 



(5.9) 



holds. This condition follows from the fact that, the char- 
acteristic area enclosed by the relevant trajectory should 
not exceed X'j^. If Eq. ( |5.9| ) is not fulfilled, the trajectory 
should turn back at the distances much smaller than Itr- 
The probability of such an event is determined by the 
optimal configurations consisting of a small number of 
scatterers andp— thus, separation of the diffusion region 
is not possibleEy. In all the sub sequent calculations, we 
assume that the condition ( ^.9[) is met. 

In the diffusion region, the Cooperon satisfies the equa- 
tion 



-iuj- DlV 



(Cgi)=5(p). (5.10) 



At a ^ p ^ Itr, the Cooperon Cgi ceases to depend on 
p, and we have with the logarithmic accuracy 



(Cgi) 



1 



AttD 



In 



UJTtr 



-Y 



D 



(5.11) 



where dimensionless function Y (x) is given by0 

y (x) = * Q + +ln4a;, (5.12) 

and ^'(x) is the digamma function. 

The solutio n in t he Lyapunov region with the bound- 
ary condition (5.11) can be represented in a form similar 
to Eq. (pTO| ) 



(Cgi) 



Wc(^;0, p) 
AttD 



In 



LUTtr 



Y 



D 



(5.13) 



Here, function Wc is related to Wc by Eq. (4.11), how- 
ever, the equation for the latter function, see Eq. ( 4.33 ) 
is modified: 



d_ 

di 









VF ^ 


'^dp 





d dU d 
aR ~ 9R ' 9P 



(5.14) 



PFdR\ ' 



Wc = 0,. 



Equa. tion (5.14) is supplied with the boundary condition 
( Fl3|) . 

Now, we will show that this modification does not af- 
fect functi on W g in the Lyapunov region provided that 
condition (5.!;) holds. Thus, the renormalization func- 
tion T{lo) is not affected by the magnetic field. In order 
to demonstrate this we use the following arguments. The 
effec t of t he extra in comparison with Eq. (|4.33| ) term in 
Eq. (5.14) can be taken into account by multiplying func- 
tion W±_ from Eq. ( ^.34| ) by the factor exp {2iA{t)/\y), 
where A{t) is the area enclosed by the trajectory in the 
Lyapunov region and it is given by 



A{t) 



dtiVF{tl)p{tl). 



(5.15) 



Let us estimate the maximal value of area A. In the Lya- 
punov region, the distance p does not exceed the char- 
acteristic scale of the potential a. In the vicinity of the 
boundary of the Lyapunov region p depends exponen- 
tially on time p{t) ~ ae^* (here time t < is counted from 
the moment of arrival of the trajectory to the boundary 
of th e Lya punov region). Substituting this estimate into 
Eq. ( 5.15 ), we obtain 



maxj^l ~ avp/X ^ Ifr- 



(5.16) 



Comparing estimate ( 5.16 ) with the condition (5^), we 
conclude that |^| Xjj and, therefore, the magnetic 
field has no effect in the Lyapunov region. 

Thus, final formula for the weak localization correction 
in the magnetic field H reads 



A(7{H,Lu) - Acr(0,w) 



r{Lu)Y 



D 



(5.17) 



where functions F and Y are defined by Eqs. (1.2) and 
( ^.12| ) respectively. It is worth noticing that the effects of 
the phase relaxation, see Eq. (5.4), and of the magnetic 
field on the renormalization function are different. This 
is because the effect of the phase relaxation is determined 
by the time the particle spends in the Lyapunov region, 
which is significantly larger than rt^, whereas the effect 
of the magnetic field is governed by the area enclosed by 
the trajectory in the Lyapunov region which is always 
much smaller than Ij,.. 

For the weak magnetic fields, Xjj ^ D/niax{u!,T^ 
we obtain from Eq. 

Aa{H) -Acr(O) = 



(5.17) 



F (cj) F cj + 



1/2 



(5.18) 
n 2 



(1 - ILUT^) Xjj 



The study of the frequency dependence or temperature 
(via Tip) of the magnetoresistance may provide an addi- 
tional tool for measuring the Lyapunov exponent. 
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VI. WEAK LOCALIZATION IN THE BALLISTIC 
CAVITIES 

In this section we study how the Lyapunov region af- 
fects the weak locahzation correction in the balhstic cav- 
ities. At zero- frequency and T^p —t oo this problem was 
studied in Refs. [|-|lO). 

For the sake of simphcity, we restrict ourselves to the 
case of zero magnetic field H = and concentrate upon 
the dependence of the weak localization correction to the 
conductance Ag of a ballistic cavity on frequency w and 
phase relaxation time r^. The effect of the magnetic field 
on the weak localization was studied in Ref. iM]. 




FIG. 4. Schematic view of the ballistic cavity "B" between 
two reservoirs "L" and "R". 

Let us consider the system consisting of three cavities, 
see Fig. |[ connected by channels. The size of the central 
cavity ( "B" in Fig. ^ is much smaller than that of the 
outer cavities ( "L" and "R" in Fig. ^ which act as reser- 
voirs. The conductance of the system is controlled by the 
channels so that their widths dL,R are much smaller than 
the characteristic size of the central cavity, dL.n. o-- 
We assume that the motion of an electron in the channel 
still can be described by the classical Liouville equation, 
which implies (Il^r ^ A^. 

Because of the inequality d^^R <^ a, the time it takes 
to establish the equilibrium distribution function in the 
cavity is much smaller than the escape time. (The equi- 
libration time is of the order of the flying time of the 
electron across the cavity.) Under such conditions the 
classical escape times from the cavity through the left 
(right) channel T/^(fl) are given by 



tl{r) 



1 



'-^Join-di) VLiR)n-di^(^y- 



dn 
2 



'■L{R)Vl{R) 



2Ab 



(6.1) 



where Ab is the area of the cavity, the linear integration 
is performed along the narrowest cross-section of the cor- 



responding channel, dl^^^^ is directed outside the cavity 
"B" normal to the integration line, and vl( r) ^re the 
Fermi velocities in the contacts. Equation (6.1) corre- 
sponds to the classical Sharvin formulacll for 2D case and 
the escape times are related to the classical conductance 
of a single channel gL{R) by 



9l{r) 



tl(r) 



(6.2) 



If the external bias eV{t) is applied to, say, left reser- 
voir, (the right reservoir is maintained at zero bias), the 
electric current / from the left to the right reservoir ap- 
pears. This current is linear in the applied bias: 



I{t) = ~QL{t) = / dt'g{t-t')V{t'), 



9{t) 



duj 
— f 
27r 



(6.3) 



where Qj^ is the charge of the left reservoir. Relation 
Eq. ( |6?^ ) defines the conductance of the system g{uj). 
Performing actual calculations in Eq. ( |6.3| ), one has to 
take into account the condition of the electroneutrality 
in the cavity "B", Qb = 0. The electroneutrality re- 
quirement is valid at times larger the characteristic time 
of the charge relaxation. This time Tc can be estimated 
as Tc — CB/(max gL.R), where Cb is the capacitance of 
the cavity. Using estimate Cb ~ a and formulas (|6.l|), 
(F^ i we find Tc ~ T//aB/(max dL,R), where r// = a/vp 
is the flying time of the electron across the cavity, and qb 
is the screening radius in 2D electron systems. For wide 
channels dL.R^ ^ Xp — clb, we have Tc ^ r/;. We are 
interested in the dynamics of the system at time much 
larger than the flying time and, therefore, we can assume 
that the electroneutrality holds. 

Then, the standard linear response calculations enable 
us to relate the conductance g to the diffuson T> defined 
in Sec. |l[ The charge response in ith cavity, Qi to the 
applied biases V{t), VB{t) = V, Vee"'"* can be expressed 
by means of the polarization operator as 

Q, = e^jdvidv2 n (l^ ; r 1 , r2 )0. (r 1 ) (r2 ) + Us 0i3 (r2 )] , 

(6.4) 

where function Oi (R) equals to unity if vector R lies in 
the ith region (i — L, i?, B) and equals to zero otherwise. 
The potential Vb is to be found self-consistently from 
the electr oneu trality requirement. Subs titut ing Eq . (|3.4| ) 
into Eq. (6^) and making use of Eqs. (3.5), and ( p.l2[ ), 
we obtain with the help of definition (|6.3|) 

Vn 

Vll{uj)+Vlb{uj)^ 



g{ui) — se V <y^iujAL + ^ 

dnidn2 dRi (iR2 



&,(Ri)0,(R2) 



xI?,y(a;;ni,Ri;n2,R2), 



(6.5) 
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where Ai is the area of the corresponding region (z = 
L,R,B). 

The electroncutrahty condition, Qs = 0, gives us 
the equa tion for the potential of the cavity Vb- Usin g 
Eq. (jej) for i B, wc find with the help of Eqs. (Q, 
(O and (|3.12D 



iuVLBito)V + [Ab + iLoVBBii^)] Vb - 0. (6.6) 



A. Classical conductance 



Let us first calculate the classical conductance gd of 
the system. We consider the frequencies u, much smaller 
than the flying time of the electron in a cavity. Assuming 
that the motion in the cavity is ergodic and the areas of 
the reservoirs are large, ^l(_r)/^_b ^ ^'''l(r)-, we obtain 
that the diffuson changes only within the channels. For 



"Dij from Eq. (6.5) we find 



Ab 



1 1 1 



-iijj + — tb tl tr 



T)" (f A — J -X- ^ BB 

(TjW) 



J=L,R; 



J=L,R; 



(6.7a) 
(6.7b) 
(6.7c) 
(6.7d) 



Equation (6.7a) describes the exponentially decaying in 
time probability to find the electron in the c avity "B" 
if it started in this cavity. First term in Eq. ( |6.7b ) cor- 
responds to the classical correlator of the jth reservoir 
disconnected from the cavity, the second term describes 
the finite probability for the electron to enter cavity "B" 
from jth reservoir, and the third term corresponds to 
the process in which an electron from jth reservoir visits 
the cavity once and then comes back. Equation (6.7c) 
gives the probability for the electron to appear in the jth 
reservoir starting from the cavity. Finally, Eq. ( ^.7dD is 
the probability for an electron to get from the left to the 
right reservoir. 



Substituting Eqs. (6.7) into Eq. (6.6), we find that 
the bias of the cavity Vb does not depend on frequ ency , 
Vb = gbl [gh + fffl)- Then, by substitution Eqs. (6.7) 
into Eqs. (6.5), we obtain with the help of Eq. (|6.2|) 



9ci 



9LgR 
9L + 



(6.8) 



in agreement with the Kirchhoff law. ft is worth men- 
tioning that the result ( |6.8[ ) at w = can be obtained 
without requirement of the electroneutrality. 



B. Weak localization correction 

In order to calculate the weak localization correction 
to the conductance /S.g{uj) we have to find the c orre ction 
t o th e classical correlator AP and then use Eq. ( |6^ ) and 
( |6.6|). For such calculation, it is most convenient to use 
Eq. ( [3.28 ) . Our strategy will be analogous to the one we 
used in Sec. |^ for the calculation of the correction to 
the conductivity. 

Integrating both sides of Eq. ( 3.28|) over the coordi- 
nates 1 and 2 within the regions specified by ^?-functions 
in Eq. (6.5) and using obvious relation I?(l, 2) = I'(2, 1), 
we obtain 



AP 



-,(2) 



J' 
dl 



ji ) 



(6.9) 



I?O(l)0,(Ri)+2??(l)0,(Ri) 



c°(i,i) 



AV^' = / dll?,"(l)l?°(l) 



2iu! — Li 



1 

Tq d(j)l 



where 



I?°(l)= / d20,(R2)I?°(2,l) 



C°(l,l) 
27rv '' 



(6.10) 



Here, we use the short hand notation I = (ni,R;), inte- 
gration over the phase space on the energy shell is defined 
as dl = (in;dR;/27r, and the time reversed coordinate / is 
given by [= (-ni,R;). 

It is noteworthy, that the quantum correction AVij 
satisfy the charge conservation condition 



E 

j=L,R,B 



AP,;, 



0, 



-L, i?, B 



(6.11) 



which can be easily proven with the help of the relation 
^.X'O(l) = -J- and Eq. (|^). Equation ( loTl ) en- 
ables us to consider only non-diagonal elements of AVij 
which is technically easier. 



Analogous to the discussion in Sec.lV, we assume that 
the Cooperon part of the expression can be calculated in- 
dependently of the diffuson part. This is because the clas- 
sical trajectories corresponding to these quantities tra- 
verse essentially the different regions of the phase space. 

First, we use this assumption to evaluate contribution 
from Eq. (6.9). We notice that the classical tra- 
jectory can close only inside the cavity. Therefore, the 
Cooperon C(l, 1) also exists only inside the cavity "B". 
For the calculation of the diffuson, we notice that at times 
much larger than the flying time across the cavity Tfi the 
position of the electron and its momentum is randomized. 
It suggests to use the approximation 



« ^ / dlP°(l)0B(Ri) = 

^B J ^B 



(6.12) 



if vector Ri lie s in side the cavit y. H ere, function is 



defined in Eq. (|6.7D . Using Eq 



y. M e 



we obtain 



B"iB 



(6.13) 
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where the average inside the cavity is defined as 

{...) = j- J dWB{Ri) .. . . 
Let us turn to the calculation of the contribution AP, 



(2) 

As we already saw in Sees. || and EV, two diffusons can 
not be averaged independently, because motion of their 
ends are governed by the same potential during period 
tE 3> Ttr- On the other hand, the randomization of the 
motion of the center of mass occurs during in a time in- 
terval of the order Tfi. Therefore, we can approximate 



P,°(1)I?"(1)«0b(Ri)(I?,°(1)I??(1)). 



(6.14) 



Expression (6.14) is written in the lowest order in small 
parameter Ab / ^L.R^ and we exclude from our consid- 



eration cases i 



3 



L, R. In the latter cases, there 



are also non- vanishing contributions in Eq. ( 6.14 ) cor- 
responding to the coordinate Ri in the reservoirs L or 
R. This would require more careful investigation of the 
behavior of the diffuson in the channels. We, however. 



simply bypass this difficulty by utilizing identity (6.11) 
for the calculation of the diagonal elements AVll and 
AVjfii. Using the approximation (6.14), we find 



ADg' = (I>,»(l)r;(l))/dl9B(R, 



2iLU — Li 



C°(l,l) 
(6.15) 



Liouvillcan operator Li in the second term of the RHS 
of Eq. ( lel^ ) is the total derivative along a classical tra- 
jectory and, therefore, it can be reduced to the linear 
integrals across the channels 



dl ^^B(Ri)LiC°(l,i) 



dni 
"2^ 



ni 



d£^ 



ni 



(6.16) 



fFC°(l,l) 



wher e th e linear integration is defined similar to that in 
Eq. (|0|). Then, we notice that a classical trajectory can 
close only inside the cavity. Therefore, only the Cooper- 
ons with the initial momentum directed inside the cav- 
ity exist. Assuming that the randomization of the mo- 
mentum direction occurs only inside the cavity and con- 
sidering the times much larger than the flying time, we 
conclude that the Cooperon in the contact vanishes if its 
moment rii directed inside the cavity and for the moment 
directed outside the cavity Cooperon coincides with its 
value inside the cavity, C(l, 1) = 6* (rii • (i^L(fl)) (C(l, 1)) 
for the coordinate Ri located in the left or r ight channel 
respectively. This enables us to reduce Eq. ( |6.16D to the 
simple form 



dl 9b(Ri)LiC"{1,1) 



TB 



■(C(1,T)), 



.17) 



where the tot al esc ape time tb is defined in Eq. (3.7a). 
Deriving Eq. ( 3.17| ), we use the definition of the escape 



times (^^). Arguments above are essentially equiva- 
lent to those-jn the derivation of the classical Sharvin 

conductanceEll. 

Combining formulas ( 6.13 ) , ( 6.15| ) , ( 6.17 ) and ( |6.9|) , we 
obtain 



2711^ 



(6.18) 



2^u;--)AB{Vnl)'D]il)) 
Tb J _ 



We reiterate that Eq. (6.1?) is not applicable for the case 
of i = j = L, R. In order to find the diagonal elements 



AVll and AVrr, one has to use the identity (|6.11 ). 

The calculation of the corresponding averages 
(C°(l,i)) and (I?,"(l)X'°(l)_)Js performed along the lines 



of the derivations in Sec. [V. In the calculation of the 



Cooperon the only change is in the expression (4.3) for 
the Cooperon outside the Lyapunov region 



C{4>,p) 



1 



Ai 



1 



(6.19) 



which is analogous to Eq. (3.7a) The solution for the 
Cooperon in the L yapun ov r egion is analogous to one 
presented in Sees. IV A and tVE. The calculation of 



function w may be performed for the cavity disconnected 
from the reservoirs, provided that the condition Ie tb 
holds. As the result we obtain 



(C(i,i)) = 



w(w,0,0) 



Ai 



1 

TB 



(6.20) 



In the calculation of the product of two di ffuson s 
(X'°(1)X'°(1)) the change should be made in Eq. ( |4.25D . 
The reason for thi s is t hat the integration over R, n for 
the reducing Eq. (4.22) to Eq. (4.24) is performed now 
only inside the cavity. As a result, one more term 



■0s(R)LcMi(l,2;n,R;,/),p) « - ' '^'^^ 



27r 



TB 



[cf . with the deri vatio n of Eq. ( 6.17])] ha s to be added to 
the LHS of Eq. ( [4.241 ). Equation ( [l.25D , then, acquires 
the form 



Mi(l,2;p, 



(I?°(1;2)) + (I?°(1;2)) 
-2iuj + 



and we obtain instead of Eq. (4.28) 



Ab{V^{1)VUi))^w{lo; 0, 0)^'^^ 



Ab 



+ 



1 - u;(w;0,0) 
-2iuj + 

TB 



[D%5,B+V%5,B^. (6.21) 
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where f unctio ns V^g are given by Eqs. (3.7a) and (6.7c). 
Result ( 6.21 ) is no t applicable for i = j = L,R cases. 
Deriving Eq. ( 3.21 ), we used Eq. ( 6.12| ) for the average 
of the single diffuson (2?° (1)) 

Substituting Eqs. mM) an d (|6.2lD into Eq. ( |6.18| ), we 



obtain with the help of Eqs. ( |6.7D and ( [4.32|) 



AVbb H = 



1 

Tb 



2'KV 



1 

Tb 



AVjB = ^'Dbj = - — AVbb, 



AV 



1 — 2iujT 



LR. 



J=L,R; 



AVbb H . 



(6.22a) 

(6.22b) 
(6.22c) 



Corr ectio ns AVjj for j = L,R are found with help of 
Eq. ( |6.11 ) and they are given by 



2iLUTB — 1 



— ] AVbb- 

Ti 



(6.22d) 



Substituting Eqs. ( |6.22a| ) and ( |6.22c| ) into Eq. (|6^, 
we observe that the voltage in the cavity Vb does not 
acquire any quantum corrections, Vb — VgiJ {9r + 9l)- 
Finally, substituting Eqs. ( |6.22 ) into Eq. ( |6.5| ) and using 
Eq. (3.2), we obtain the final result for the frequency de- 
pendent weak localization correction to the conductance 
of the ballistic cavity 



Ag{u) 



9 LOR 



V{lo) 



1 — iujTB 



(6.23) 



where the total escape time tb is defined in Eq. (6.7a) 
We emphasize that Eq. ( ^.2S ) at zero frequency can be 
obtained without electroneutrality requirement. 

Equation ( |6.23| ) is the main result of this section. At 
= 0, this result agrees with the findings of Ref. [ ||] 
in the limit of large number of quantum channels in the 
contact. We are aware of neither any calculation at finite 
frequency nor of a description of the role of the Ehren- 
fest time in the conductance of the ballistic cavities. The 



renormalization function r(a;) in Eqs. ( 6.23 ) describes 
the effect of the Lyapunov region o n th e weak localization 
and it is given by Eqs. ( |l.2| ) and (L3). Analytic calcula- 
tion of the Lyapunov exponents A, A2 ~ t7^ for the ballis- 
tic cavity is a separate problem and it will not be done in 
this paper. It is assumed in Eq. ( ^.2S ), that the condition 
ts tb holds. The result for the opposite limit, (which 
corresponds to the exponentially small Planck constant), 
is o btain ed by substitution r(a;) Tlu -\- i/tb) in 
Eq. ( |6.23| ) and the weak localization correction turns out 
to be suppressed by the factor exp {—2tE/TB). 

The finite phase relaxation time is take n int o ac- 
count by substitution u; uj + i/T^p in Eq. (6.2C). At 



Tip S> i_E, the result for dc- conductance agrees with the 
result of Ref. [ |2^ . We obtained for w = 



Ag: 



se 



9L9R 



2nn 



(gL+gR) 



(6.24) 



exp 



tE 



A27 



TB 



1- 



A2 



tbj 



where Ti is the time it takes for an electron to be scattered 
inelastically or to escape the cavity, 



1 

n 



1 

TB 



1 

Ti 



Usually, the Ehrcnfest time is is much smaller than the 
escape time tb- In this case, one can immediately see 
the dramatic crossover at the temperature dependence 
(usually Tip is a power function of temperature, see Ref. [ 
[1^). If at 3> ifi, the dependence on temperature is 
a power law, with the increase of the temperature the 
change to the exponential drop occurs. Thus, the study 
of the crossover in the temperature or frequency depen- 
dence of the ballistic cavities may provide the information 
about the values and the distribution of the Lyapunov 
exponents in the cavity. 



VII. CONCLUSION 



In this paper we developed a theory for the weak lo- 
calization (WL) correction in a quantum chaotic system, 
i.e. in a system with the characteristic spatial scale of 
the static potential, a, being much larger than the Fermi 
wavelength, \p. We showed that for the quantum chaos, 
new frequency domain appears, ^ <C t, 



[tE 

is the Ehrenfest time, see Eq. (1.2)] where the classical 
dynamics is still governed by the diffusion equation, but 
the WL correction deviates from the universal law. For 
the first time, we were able to investigate frequency de- 
pendence of th e W L correction at such frequencies, see 
Eqs. (1.1) and ( |1.2D , and to find out how the fundamen- 
tal characteristic of the classical chaos appears in the 
quantum correction. At lower frequencies, lo <^ t^^, we 
proved the universality of the weak localization correc- 
tion for the disorder potential of an arbitrary strength 
and spatial size. 

These results may be experimentally checked by study- 
ing the frequency or temperature (via T,p) dependence of 
the weak localization correction {e.g. negative magne- 
toresistance). Indeed, at the low- frequency or temper- 
ature, the conventional dependence should be observed. 
This dependence is rather weak (logarithmical for large 
samples and a power law for the ballistic cavities). With 
the increase of the frequency or temperature, the depen- 
dence becomes exponential; such a crossover may be used 
to find the Ehrenfest time, tp and thus extract the value 
of the Lyapunov exponent. The parameters of the ballis- 
tic cavities studied in Ref. [[t) are a ~ 1/Ltm, Xp — 400A, 
so that ln(a/AF) — 3. We believe, however, that the size 
of the ballistic cavities may be raised up to the mean free 
path ~ 17/im; Ehrenfest time in this case would be ap- 
preciably larger than the flying time, ln(a/A_F) ~ 6, and 



20 



the characteristic frequency uj = for this case can be 
estimated asu;~5xl0^s~^. Measurements of frequency 
dependence of the WL correction in quantum disorder 
regime were performed in Ref. [ at frequency as high 
as 16.5 GHz. Thus, the measurement of the Ehrenfest 
time in the baUistic cavities does not seem to be unreal- 
istic. 

We expect that the effects associated with the Ehren- 
fest time may be found also in optics. They may be ob- 
served, e.g. in the dependence of the enhanced backscat- 
tering on frequency of the amplitude modulation to, This 
dependence should be still given by our function r(cL;) 
with Xp being replaced with the light wavelength. 

We showed that the description of the intermediate re- 
gion ^ w ^ Tj~^ can be reduced to the solution 
of the purely classical equation of motion, however, the 
averaging leading to the Boltzmann equation is not possi- 
ble because the initial and final phase cells of the relevant 
classical correlator (Cooperon) are related by the time in- 
version. Therefore, the initial and finite segments of the 
corresponding classical trajectory are strongly correlated 
and their relative motion is described by the Lyapunov 
exponent and not by the diffusion equation. We took 
this correlation into account, showed that it is described 
by the log-normal distribution function and related the 
Ehrenfest time to the parameters of this function. 

Because the description by the Boltzmann equation 
was not possible, we derived the lowest order quantum 
correction to the classical correlator in terms of the solu- 
tion of the Liouville equation, smeared by the small angle 
diffraction, see Eq. ( 3.27bD . The derivation was based on 
the equations of motion for the exact Green functions 
and did not imply averaging over the realization of the 
potential. 

Closing the paper, we would like to discuss its relation 
to the other works and to make few remarks concerning 
how the Ehrenfest time appears in the level statistics. 

First, we notice that, though quite popular in the clas- 
sical mechanics and hydrodynamics, the Lyapunov expo- 
nent very rarely enters in the expressions for observable 
quantities in the solid state physics, see Ref. [ The 
possibility to observe the intermediate frequency region 
t~^^ <^ uj <^ Ttr appeared only recently with the techno- 
logical advances in the preparing of the ballistic cavities 
and that is why the region has not been studied system- 
atically as of yet. Let us mention that the importance 
of the Ehrenfest time in the semiclassical approximation 
was noticed already in Ref. [ |lj| where it was shown that 
the method of quasicJassical trajectories in the theory 
of superconductivityllj fails to describe some non-trivial 
effect at times larger than tp which was calculated for 
the dilute scatterers. The term "Ehrenfest time" for the 
quantity ( |1.3| ) was first introduced in Ref. [ Q The rele- 
vance oHe in the tiieory of weak localization was empha- 
sized by ArgamarJlJ, however, he focused only on times 
much larger than the Ehrenfest time. 

The universality of weak localization correction at 
small frequencies was known for the case of weak quan- 



tum impurities^ and for the ballistic cavitieaa. We are 
not aware of any proof of the universality for the disor- 
der potential of arbitrary strength and spatial scale. 

The description of the quantum corrections in terms of 
the nonaveraged solutions of the Liouville equatipH was 
developed in by Muzykantskii arid, KhmelnitskiiEs and 
more recently by Andregw et. alr^, who suggested the 
effective supresymmetricEH action in the ballistic regime. 
In Ref. [ |2^, the supersymmetric action was written in 
terms of the Perron-Frobenius operator which differs 
from the first order Liouville operator by the regular- 
izator of the second order. This regularizato r is s imilar 
to the angular diffusion term, c>c l/r^ in Eqs. ( 3.27 ). The 
authors mentioned that all the physical results can be ob- 
tained if the limit of vanishing regularizator is taken in 
the very end of the calculation. Our finding indicate that 
the time it takes for the quantum correction to reach its 
universal value is oc In (r^). Thus, at any finite frequency, 
the limit — > cx) can not be taken and the regulariza- 
tor in the supersymmet ric a ction should be assigned its 
physical value, see Eq. ( 3.26 ). 

In principle, our formula for the weak localization cor- 
rection (3.27b) can be derived using the supersymmetry 
technique. However, our approach seems to be techni- 
cally easier and more physically tractable for the calcu- 
lation of the first order weak localization correction. We 
believe that the supersymmetry may serve as a powerful 
tool for the investigation of the effect of the Ehrenfest 
time on the higher order corrections and on the level 
statistics. 

It is generally accepted that the level statistics at law 
energies is described by the Wigner-Dyson distributionEa. 
For the simall disordered particle it was first proven 
by EfcteacZl and for the ballistic cavities by Andreev 
et. alx^.—. For the quantum disorder, Altshuler and 
Shklovskiita showed that the universal Wigner-Dyson 
statistics breaks down at the Thouless energy. For the 
ballistic cavities the universal statistics is believed to be 
valid up to the energies of the order of the inverse flying 
time r//, at smaller energies s the corresponding correc- 
tions are small as STfi. We, however, anticipate devia- 
tions at the parametrically smaller energies of the order 
of t^^, and the corrections of the order of stp at energies 
s < t^^ 

Let us consider for concreteness the correlator of the 
density of states 'R{s) = {p{e)p{e + s)) — {p{e))^, where 

p(e) = TrS — . For the orthogonal gaussian ensem- 
ble the random matrix theory yieldsEl TZ{s) = — (tts) + 
(l -I- cos^ tts) / (tts) * -I- . . ., where s ^ 1 is measured in 
units of mean level spacing. We expect, that the first 
term in this expression is not affected by the presence 
of the Lyapunov region, whereas the following terms are. 
In the supersymmetric approachcZI this follows from the 
fact that the first term arises from non interacting dif- 
fuson modes whereas all the others come from the inter- 
action of these modes. Such interaction is analogous to 
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the one giving rise to the weak localization which was 
shown to have the frequency dispersion des cribed by the 
renormalization function r(a;), see Eq. (p^). We believe 
that the same renormalization factor will appear in all 
the effects associated with the coupling of the diffuson- 
Cooperon modes. 
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APPENDIX A: LYAPUNOV EXPONENT FOR 
THE WEAK SCATTERERS 

We consi der explicitly the case where the potential U 
in Eq. (4.33) is weak and its distribution function is Gaus- 
sian. For the sake of simplicity we neglect the angular 
diffusion due to the quantum impurities in the Lyapunov 
region, becau se thi s diffusion does not affect values of 
A, A2, see Sec. IV E . In this case it is more conveni ent no t 
to follow the general procedure outhned in Sec. ( I V E ) , 
but to make use of the small parameter a/ltr first. Con- 
sidering the disorder potential in the second order of the 
perturbation theory, we obtain for independent on R, n 
part of the function W± : 



d_ 

dt 



Vf4> 



d_ 

dp 



2 

Ttr a<p'^ 



(Al) 



where the transport life time is given by 



1 



1 



dx{dyU{x,Q)dyU{Q,Q)), (A2) 



Ttr 4:eFPF 

and the dimcnsionless function £ is defined as 

C^^dx{dyUix,p)dyUio,o)) 



£{p) 



1 



/_°°^dx(9,C/(a;, 0)9,(7(0,0))- 



(A3) 



In Eq. (Al), we ass umed only and lifted the other 

assumption of Eq. ( 4.33| ) p <^ a. If p <C a, we expand 



£ in Taylor scries, £{p) ~ p /2a , which rigorously de- 
fines length a in this case, and we arrive to the equation 
describing the Lyapunov region for the weak disorder po- 
tential 



d d p^ 
— - vf<Pt, 2 

at op Ttrtt'^ 
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Wa^ = 0. 



(A4) 



It is worth noticing, that our approach is equivalent to 
one involving the multiplication of the vector (p, (p) by a 
Monodromy matrix after each scattering event. Equation 
( A4) is valid because each Monodromy matrix defined on 



a time of the order a/vp is close to unit matrix. Other- 
wise, the last term in the brackets in Eq. (A4) becomes 
an integral operator. 

After introduction of new variables 



= In- 



V = 



A /I \ 1/3 



t / Itr 
Ttr \ a 



2/3 



equation (A4) acquires a simple form 



d_ 
97 



d 



d 



oz oy ay^ 



W_L = 0. 



(A5) 



(A6) 



We are interested in the case when function W± 
changes slowly as the function of z. Corresponding gra- 
dient is small, and we can employ the procedure similar 
to the reducing the Boltzmann equation to the diffusion 
equation. Let us represent function W± as 

W^{t- z, y) = WI{t; z) + W\t- z, y), Wi « W±. 

(A7) 



Substituting Eq. ( |A7| ) into Eq. ([Aq), multiplying the re- 
sult by function g{y): 



ay 



dy^ 



g{y) = 0, / dyg{y) = 1 



(A8) 



and integrating over y we obtain 
d 



Ot dz 



dy [iy-P)9iy)w'iy)] =0, 

(A9) 



where the numerical coefficient /3 is given by 

P= fdy yg{y) (AlO) 



(All) 



and function W can be written as 



= h{y) 



dy dy"^ 



Shift of W^^ by an arbitrary constant does not affect 
the results, see Eqs. ( A13 ) and ( [AlOj ). Substituting 
Eq. (All) into E lq. ( A9), we obtain in accordance with 
general formula (4.47) 



OT oz az^ 



(A12) 



WI{t-z)^Q, 
where the numerical coefficient /32 is given by 

P2= [dy {y-p)g{y)h{y). (A13) 
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Comparing Eqs. ( |A12| ) and ( |A5| ) with Eq. ( [4.47D , we find 
for Lyapunov exponent A and its deviation A2 



Ttr 



2/3 



W2 fk^ 
Ttr \ a 



2/3 



(A14) 



Simple calculation of the numera l coef fic ient s /3 , /32 i s car- 
ried out with the help of Eqs. ( |A10D , (^), ( |A13D and 
( |All| ) with the final result 



e-y''/^ r dyi evf/^ 

J— 00 ^ J— 00 



0.365, 



(A15) 



J—oo^y 



/32 = 



r dy e-y'/^ r dyi evU^ 

J— 00 J —OQ 



0.705. 



In order to avoid any confusion, let us notice that the 
log-normal distribution function can not be used to find 
the averaged moments of the coordinates p, (/) and it is 
sufficient only for the calculation of the low moments of 
the logarithm of the coordinates. 
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^ This result, strictly speaking, is applicable for the case of 
the strong scatterers where the transport and elastic mean 
free paths are of the same order. For the case of weak scat- 
terers one should use vf/A instead of Itr- Here A is the 
Lyapunov exponent, see Appendix ^. 
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